Influence of the atomic scale inhomogeneity of the pair interaction on the local pair 
formation and density of states in high-T c superconductors. 
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The influence of the atomic-scale inhomogeneities of the pairing interaction on the superconduct- 
ing order parameter distribution and the LDOS is studied in the framework of mean-field BCS theory 
for two-dimensional lattice model. It is found that the ratio of the local low-temperature gap in 
differential conductance to the local temperature of vanishing the gap 2A g /T p can take large enough 
values compared to the homogeneous case. This ratio practically does not depend on the location in 
the sample and is independent on the concentration of local pair interaction perturbations in wide 
range of concentrations. The obtained results could bear a relation to the recent measurements by 
Gomes et. ah. 
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Nanoscale inhomogeneities have been widely ob- 
served in the high-temperature superconductor 
Bi2Sr2CaCu2C>8+ x (BSCCO) and have generated 
intense interes t 2 ' 3 ^ 4 '^' 6 "' 7 '. In particular, the spectral 
gap in the local density of states (LDOS) has been 
investigated by scanning tunneling microscopy (STM). 
It was found that the gap varies by a factor of 2 over 

o 

distances of 20 — 30 A. These observations have been 
made primarily in the superconducting state. Several 
scenarios have been proposed to understanding this 
electronic inhomogeneity. First of all, it was speculated 
that poorly screened electrostatic potentials of the 
dopant atoms vary a doping concentration locally, 
giving rise to the gap modulations^&i -. Alternatively, 
these inhomogeneities are associated with a competing 
order paramete r 11 : 12 : 13 ! 14 : 15 : 16 . Further, the positive 
correlations between the inhomogeneities and positions 
of the dopant atoms have been observed by STM on the 
optimally doped BSCCO 6 . After that it was proposed 
by Nunncr et al~Z- that the dopant atoms modulate 
the pairing interaction locally on the atomic scale. The 
LDOS calculated in the framework of this model is 
in good agreement with the key characteristics of the 
experimental data. In addition, finite-temperature order 
parameter evolution was investigated and the specific 
heat was calculated within the same models. 

On the other hand, it is well known that in the high- 
Tc superconductors a partial gap in the LDOS exists for 
a range of temperatures above There is no con- 

sensus up to now if this gap is due to pairing with- 
out phase coherence, a competing order or proximity to 
the Mott stat o 20 : 21 ' 22 : 23 . The inhomogeneities described 
above complicate the situation. Only very recently the 
spatially resolved STM measurements of gap formation 
in BSCCO samples were performed 1 -. For a range of dop- 
ing from 0.16 to 0.22 they have found that gaps nucle- 
ate in nanoscale regions above T c and proliferate as the 
temperature is lowered, evolving to the spatial distribu- 
tion of gap values in the superconducting state. It was 
observed experimentally that overdoped and optimally 
doped samples have identical gap-temperature scaling ra- 



tios, which together with the fact that in the overdoped 
samples pseudogap effects are believed to be weak or ab- 
sent, allowed Gomes et.al to interpret the gaps above T c 
as those associated with pairing. The most striking ex- 
perimental observation is that, despite the inhomogene- 
ity, every pairing gap develops locally at the temperature 
T p following the relation 2A/T p — 7.9 ±0.5 in wide range 
of doping from overdoped to optimally doped samples. 
So large and independent on the size of local gap val- 
ues ratio 2A/T p seems to be different from the expecta- 
tions based on the BCS theory and its strong-coupling 
extension, where it is in the range 3.5 — 5 and becomes 
dependent on A. 

In the present paper we show that if the pairing inter- 
action is modulated locally on the atomic scale, as it was 
proposed by Nunner et. aS^, the ratio 2A/T p strongly 
increases in the framework of conventional BCS theory. 
This is in sharp contrast to the influence of the potential 
disorder, which, as was demonstrated^ 4 ., can only dimin- 
ish this quantity. It is found that the ratio 2A/T p is prac- 
tically independent on the number of off-diagonal (pair- 
ing interaction) scatterers, which is proportional to the 
doping in the model by Nunner et. al 17 , in a wide range 
of concentrations. We obtained the ratio 2A/T p to be 
position-independent in case if the off-diagonal scatter- 
ers are distributed rather uniformly and equal to 5.9. For 
comparison, the ratio 2A/T p — 4.7 for the homogeneous 
sample corresponding to the lattice parameters we use (it 
is slightly higher than the ratio 2A/T C — 4.3, because T p 
is smaller than the mean-field critical temperature T c as 
it is discussed below). These results are reminiscent of 
the measurements by Gomes et. al~ except for the fact 
that the ratio 2A/T p we calculated is about thirty per- 
cent smaller than experimentally measured one. On the 
other hand, the ratio 2A g /T p is found to be very sen- 
sitive to characteristic size and height of the individual 
pairing interaction scatterer and especially to the choice 
of the lattice parameters and can be further increased 
manipulating by these quantities. It is worth noting that 
our analysis is phenomenological and the conclusions are 
independent on the underlying pairing mechanism and 
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the particular course of the local pair interaction modu- 
lations. 

Model and method.— We consider the following mean- 
field Hamiltonian on a square lattice 

{ij),a i,<r (ij) 

(1) 

where Ci a (cl a ) stands for an electron annihilation (cre- 
ation) operator at site i with spin a. 53 to) indicates 
summation over neighboring sites, tij is the hopping in- 
tegral between sites i and j. We set ty to be t = 1 
for the nearest-neighbor hopping, t' = —0.3 for the 
next nearest-neighbor hopping and t" = 0.1 for the 
next-next one. /x is adjusted to be —1 to model the 
Fermi surface of BSCCO near optimal doping. Below 
all the energies are measured in units of t. The nearest- 
neighbor d-wave order parameter (OP) should be de- 
termined self-consistently: Ay = —gij(cil c j1 — CjiCi^). 
Following Rcf. |17| we model the inhomogeneous pair- 
ing interaction by = cjb + 5g(fi + /i)/2, where gb 
corresponds to an average background interaction and 
fi = J2 S ex P( — r is/A)/ri S , where ?-j s is the distance be- 
tween the site i and the source of the pairing interaction 
perturbation. The mean-field BCS treatment is gener- 
ally believed to be appropriate for overdoped (and, to 
some extent, optimal doped) samples. So, we do not 
consider underdoped regime, where the proximity to the 
Mott state should be taken into account. 

In order to analyze the inhomogeneous pairing correla- 
tions in the framework of the mean-field hamiltonian (TTJ) 
we exploit the fully self-consistent T-matrix technique 
for Gor'kov Green's functions. The full Green's function 
takes the form 

Gij = &ij + ^ ] G ik T km G m j . (2) 
Here f km = -^(M _1 )bK m , M km = S km + 

n 

J2Gl n V nm . In this paper we only focus on the off- 

n 

diagonal self-energy inhomogeneity and therefore V km — 
(5Afe m i(T2if2. All Green's functions and T-matrices are 
4x4 matrices in the direct product of spin and particle- 
hole spaces, what indicated by the symbol ". fi and &i 
are Pauli matrices in particle-hole and spin spaces respec- 
tively. The summation is taken over all the sites, where 
the OP Afc m = A° m + SAkm differs from the background 
value A Q km . A° m is assumed to be of d-wave type, that 
is A?- 1 - = — A - = A , a and b are basis vectors of the 

Lima ii±b 

square lattice. We set the lattice constant a to be equal 
to unity. 

Local order parameter.— OP is to be calculated from 
the self-consistency equation 

A»j = gijTy^Tzj \t-ia y Gij{e n )\ , (3) 



where f_ = (f x — if y )/2. Eqs. ([3]) and ([2]) allow us to 
find the OP Ay numerically. 

We begin by considering a single perturbation of the 
pairing interaction of the form fi = exp(— r/X)/r, where 

r = i x + i 2 + z 2 . However, the results do not depend 

qualitatively on the particular form of the perturbation 
and are only controlled by its effective width and height. 




FIG. 1: (a) Dependence of the site-averaged OP maximal 
value A max on temperature for z = 0.5 and A = 0.5. g max = 
go + 8g max = 0.857, 1.022 and 1.187 from the bottom curve 
to the top one. (b) Evolution of A maa: (T) dependence with 
increasing of the perturbation size for constant g max — 1.022. 
z = 0.5, A = 0.5 (solid line); z = 2, A = 2 (dashed line) and 
z = 3, A = 3 (dotted line). Bottom and top thin solid lines 
represent homogeneous bulk dependence Ao(T) for gb = 0.692 
and gb = 1.022, respectively. 

The dependence of the site-averaged OP maximal 
value A ma * = (A , a + A ,_ a + |A og | + |A Q) _ s |)/4 on 
temperature is presented in Fig[TJ The left panel demon- 
strates the most spikier perturbation we consider: z — 
0.5 and A = 0.5. This perturbation has considerably non- 
zero value only at four bonds emanating from the cen- 
tral site. Here and below all the distances are measured 
in units of lattice constant. The different curves corre- 
spond to different heights of the perturbation: g max = 
g h + Sg max = 0.857, 1.022 and 1.187 from the bottom 
curve to the top one. The background pairing interaction 
is taken to be gb = 0.692. The thin solid line represents 
the temperature behavior of the homogeneous order pa- 
rameter without a perturbation added. We denote the 
critical temperature of the inhomogeneous sample by T x 
in order to distinguish it from the background critical 
temperature Tb- If the perturbation size is very small (a 
few bonds), as in Fig. [TJa), the corresponding T x prac- 
tically does not differ from Tb- The physical reason for 
it is that small as compared to superconducting coher- 
ence length perturbations at the mean-field level cannot 
maintain superconductivity by themselves and only do 
this due to the superconductivity in the bulk. Then the 
bulk OP vanishes, the pairing correlations in the small 
area go to zero abruptly. On the other hand, the value of 
the zero-temperature OP A max strongly enhances when 
the height of the perturbation grows. Therefore, local 
ratios /T x can considerably exceed the homogeneous 
bulk value. For the parameters we consider in Fig. []Ja) 
A /T fc = 0.58 and A max /T x = 0.76, 0.97 and 1.19 from 
bottom to top. 

When the size of the pairing interaction perturba- 



tion increases, T x grows and, correspondingly, the ratio 
Aij/T x starts to decline. Naturally, when the size of 
the enhanced pairing area becomes of the order of a few 
superconducting coherence lengths £ s , the bulk value of 
Ajj/T x corresponding to g = gj, + 5g should be restored 
at the center of the cluster. The evolution of A max (T) 
dependence with increasing of the perturbation size is il- 
lustrated in Fig. QTb) for constant g max = 1.022. For 
comparison the curve Aq(T) corresponding to this value 
of the pairing interaction is depicted in the same panel 
by the thin line. Therefore, large enough ratios A max /T x 
could be only obtained when the characteristic size of 
the OP enhanced area is considerably less than super- 
conducting coherence length. However, the anomalous 
Green's function, which enters the sclf-consistensy equa- 
tion ([3]) always has the characteristic length of the order 
of £ s ~ t/ Ao regardless of the particular parameter of the 
hamiltonian coursing the inhomogeneity. Consequently, 
the OP spacial profile follows the spacial profile of the 
coupling constant according to Eq. ([3]) except for very 
low tails extended to the distance of the order £ s due to 
the inhomogeneities of the anomalous Green's function. 
Therefore, we believe that the large enough ratio Aij/T x 
can only be obtained (at least in the framework of mean- 
field BCS theory with temperature-independent coupling 
constant) by assuming the atomic-scale inhomogeneity of 
the pairing interaction strength. 

Now we turn to discussion of more realistic situation, 
when many pair interaction scatterers are present in the 
sample. It is worth noting that in our mean-field treat- 
ment we neglect OP phase, which is a strongly fluctuating 
quantity in short coherence length high-T c superconduc- 
tors and, especially, in the inhomogeneous situation. In 
the present paper we only study local pairing and do not 
concern the global transition temperature T c , which is 
controlled by the phase fluctuations. 

We have considered a 21 x 21 sites square as a pertur- 
bation described by the T-matrix. The coupling constant 
outside this square is set to be gb = 0.692. Fig. [2] shows 
the particular example. The pairing interaction scatter- 
ers corresponding to A = 1.5 and 2 = 1 are distributed in 
the square with the concentration n = 0.078. The back- 
ground coupling constant in the square and the height of 
the individual perturbation are chosen to give the same 
average as outside the square: go = 0.45, 8g is randomly 
distributed in the range 0.70 — 1.35. 

The resulting distribution of zero-temperature OP is 
presented in Fig. [Ha). As it was described earlier for the 
single perturbation, the spacial profile of the OP inho- 
mogeneity mainly follows that one of the coupling con- 
stant. The dependence of the site-averaged Aj(T) for the 
marked sites is plotted in Fig.[^b). 

Local density of states.— STM technique measures the 
local differential conductance, and the thermally smeared 
LDOS, described by the expression 

oo 

dI/dV= f de{df{e + V)/dV) Pi {e,T) , (4) 
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FIG. 2: (Color online) (a) Zero-temperature distribution of 
the site-averaged OP. (b) The temperature dependence of the 
site-averaged OP for the sites marked by the arrows in (a). 
The bold circles show T p correspondinding to these locations, 
(c) Low-temperature gap map. (d) T v — map. 



can be extracted from these measurements. Here 
Pi(e,T) = -(l/7r)ImG£(e,T) is the local density of 
states, /(e) is Fermi distribution function and V is volt- 
age applied between the STM tip and the sample. For 
the homogeneous situation in the framework of mean- 
field weak-coupling BCS theory the distance between 
the coherence peaks 2A g in the conductance spectra 
directly connected to the OP by the simple relation 
A g ps 3.7A(T = 0) for the particular lattice parameters 
we have chosen. 

For an inhomogeneous system, where the characteristic 
size of the patch < £, there is no any direct simple rela- 
tionship between the local OP and the local gap A g . One 
can only conclude from the numerical calculations^^ 
that for regions, where the OP is enhanced from the 
background, the gap gets wider and the peak height is 
suppressed compared to the average value. It is worth 
noting that, unlike the homogeneous situation, this peak 
does not represent the maximal superconducting gap on 
the Fermi surface, but rather originates from the spectral 
weight transfer from the nearby van Hove singularity due 
to the Andreev scattering processes. Otherwise, if the OP 
in a cluster is less than that one in the background, the 
narrow and high Andreev resonant peaks develop in the 
cluster region resulting in diminishing of the gap region. 
For this reason we investigate not only the OP distribu- 
tion, but also the experimentally measurable thermally 
smeared LDOS, which has a maximum at V = A g . The 
corresponding low-temperature gap map is represented 
in Fig. [2jc) for the model sample considered above. 

Experimentally^- the temperature T p (i) of the gap dis- 
appearing for the particular location in the sample has 
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been determined using the criterion dI/dV(V — 0) > 
dl/dV (for all V > 0). Using the above criterion we cal- 
culated the distribution T p (i) from the thermally smeared 
LDOS curves. The corresponding T p -map is represented 
in Fig. [UJd). In addition, for the locations, marked in 
Fig. Et^a), appropriate T p is shown by bold circles in 
Fig. HJb) . Although in our model the OP vanishes at the 
temperature T x , which is the same for the entire sample, 
it is seen that the temperature T p is strongly position- 
dependent and considerably lower than T x for the most 
part of the sample. There are two main physical rea- 
sons for this effect: (i) LDOS is essentially nonlocal (on 
the atomic scale) quantity with the characteristic size of 
order of £ s . This fact leads to the partial averaging of es- 
sentially different gaps over the region ~ and spacial 
redistribution of the spectral weight due to Andreev scat- 
tering processes, (ii) thermal broadening of the LDOS 
further washes out the conductance peaks. It is worth 
noting that the LDOS p(e, T) taking without thermal 
smearing results in higher values of the local gap vanish- 
ing temperature. This is demonstrated in the right panel 
of Fig. 1 
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FIG. 3: (a) The probability distribution to find the particular 
value for 2A g (i)/T p (i). The bottom curve represents the same 
model sample as in Fig. [2] The top one corresponds to an- 
other example with n = 0.052, A = 2, z = 0.5, go = 0.35 and 
Sg = 0.8. The offset is for clarity, (b) LDOS (black curve) in 
comparison with the thermally smeared LDOS (gray curve) 
for a particular location and T = 0.18. 

The probability distribution to find the particular 
value for 2A g (i)/T p (i) in the model sample considered 
above is plotted in the left panel of Fig. [3] (bottom curve). 
It is seen that the distribution is quite narrow, that is this 
ratio is practically independent on the particular location 
in the sample. In the framework of our model this result 
can be only obtained if the scatterers are distributed in 
the sample randomly, but the additional restriction is 
imposed: they influence each other in the sense that the 
distance between them cannot be too small. If this re- 



striction is not imposed and the scatterers are distributed 
quite randomly, the ratio 2A g (i)/T p (i) is obtained to be 
more position-dependent. This is demonstrated for an- 
other model sample in the left panel of Fig. [3] (top curve). 
This sample differs from the considered above by the 
shape of an individual scatterer and by their concentra- 
tion. However, we have investigated a number of model 
samples and checked that the 2A g (i)/T p (i) distribution 
shape is only determined by the extent of uniformity of 
the scatterers. 

The averaged ratio 2A g /T p we calculated is equal to 
5.9. It is larger than that one for a homogeneous situa- 
tion corresponding to the same normal state hamiltonian 
2A g /T p — 4.7 but is smaller than experimentally mea- 
sured onei. However, it is obvious from the above anal- 
ysis that the underlying ratio of the OP to the critical 
temperature Aij/T x gets larger as the height of the per- 
turbation increases and its width diminishes. The man- 
ifestation of the effect in LDOS, and, consequently, the 
behavior of the experimentally measured ratio 2A g /T p 
is not so straightforward^ and strongly depends on the 
particular choice of the lattice parameters. The point 
is that for the inhomogeneous situation the peak in the 
LDOS results from the Andreev scattering processes be- 
tween the background superconducting coherence peak 
and nearby van Hove singularity, as it was already men- 
tioned above. So, the position A g of the LDOS peak 
strongly depends on the energy distance between them 
in comparison to the strength of the individual scatterer. 

In conclusion, we have studied the influence of the 
atomic-scale inhomogeneities of the pairing interaction 
strength on the OP distribution and the LDOS in the 
framework of mean-field BCS theory. It is found that 
the ratio of the local low-temperature gap in differential 
conductance spectra to the local temperature of vanish- 
ing the gap 2A g /T p can take large enough values com- 
pared to the homogeneous one. This ratio is position- 
independent in case if the off-diagonal scatterers are dis- 
tributed rather uniformly and practically independent on 
their concentration in wide range of the concentrations. 
On the other hand, it is quite sensitive to characteris- 
tic size and height of the individual pairing interaction 
scatterer and to the choice of the lattice parameters. 
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